TITLE: 
TYPE: 
ACCURACY: 
DESCRIPTION: 



METHOD OF USE: 



UNIVERSITY OF ILLINOIS 
DIGITAL COMPUTER LABORATORY 

LIBRARY ROUTIN E M 28 - 270 

Matrix Multiplication with or without rescaling (SADOl) 

Entire program 

11 decimal places (see Note l) 

The routine will form the product of two matrices, 

A., B.. = C. n . 
ij jk ik 

Features of the routine are: 

1. Matrix C can either he punched on tape or stored on 
the drum. 

2 . Matrix B can be rescaled by any power of ten prior 
to the multiplication; this enables the user to scale 
down the elements of C to prevent overflow or otherwise 
to adjust the scaling of C. 

3. A set of products of various sizes, A^'B = C^ p , 
can be produced without rereading parameters, and these 
products can all be stored on the drum for use in sub- 
sequent operations. 

k. Multiple products such as A^V 2 V ' B = C can 

be formed without punching intermediate results on tape. 

CASE I, Ay B. k - C ik St0£3 

1. Master tape 3402L 

2. Parameter tape, type T 2ko60 

3 . Data tape B 2^0J7 
k. Data tape A (with J term, sym.) 2402L 

To begin a new problem at stop 2^02L reading different 

parameters, raise the black switch. 

Matrix C is either punched on tape by rows, each element 

to "d" decimal places, or else, starting at location "D" 

on the drum, matrix C is stored by rows with a row sum 

check at the end of each row. This option is specified 

by the parameter tape. 

If matrix C is punched on tape, i, the rows of C, is not 

limited. The capacity of the computer for j and k can be 
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determined from the chart shown in the section on STORAGE 
OF MATRIX B. See also section on STORAGE OF MATRIX C. 
For alternate methods of use, see sections on CASE II, 
CASE III, and CASE IV. 

PREPARATION OF DATA TAPES § 

The elements of both B and A matrices are punched by rows 
as signed fractions with an N at the end of each row. If 
any row terminating symbol is an F instead of an N, the 
computer will stop. By raising the black switch, the read- 
ing of the matrix is resumed. (See Note 2) 
In addition to row terminating symbols, each matrix has a 
matrix terminating symbol. After the final row symbol, a 
J or an L is punched. No distinction is made between a 
J or an L at the end of a B matrix. For each,the computer 
will stop» When the black switch is raised, the computer 
will begin to read an A matrix. 

If a J terminates an A tape, when reading is resumed by 
raising the black switch, the computer will interpret the 
next tape as a parameter tape. If an L terminates an A 
tape, the previous parameters are retained, and the reading 
of a parameter tape is bypassed, 

METHOD OF USES CASE II, A^ B^ = CT P ' Stops 

1. Master tape ^hOZL 

2. Parameter tape, type T 2*4)60 

3. Data tape B^ 1 ' 2^0 J7 
k. Data tape a' 1 ' (with L term, sym.) 2ko60 

5. Data tape b' 2 ' 240J7 

6. Data tape A^ ^ (with L term, sym.) 2kO60 
etc . 

If the several C matrices are punched on tape, 50 single 

hole delays, 3 carriage returns, and a figure shift will 

separate matrices,, 

If the C matrices are stored on the drum, they will be 

stored one following the other in successive locations 

by rows with a row sum check at the end of each row. 
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PARAMETBR TAPE PREPARATION? 

Except for CASE IV, there are three parameters. The order 
in -which they must appear is as follows? 

d space number of decimal places or drum address 
+ S scaling exponent 

X directive 
Print or drum store parameter $ d. 
The parameter, d, must be followed by a fifth-hole 
character. If d < 13, the computer will punch the elements 
of matrix C to d decimal places. If d > 13, d will be 
interpreted as a drum address at which to store matrix C. 
Hereafter a drum address will be designated by a capital 
letter, D. See section called STORAGE OF MATRIX G. 

Scaling exponent, i S . 

The parameter, + S, may be omitted if no rescaling on 
matrix B is desired Otherwise, + S or -S is an integer 
exponent of 10 « Each element of the B matrix will be 
multiplied by this power of 10 before B is stored on the 
drum. This has the effect of rescaling the product 
matrix C« There are no restraints on the value of + S 
other than those imposed by the data« Of course, if the 
scaler is -12, matrix B becomes a null matrix. Similarly, 
the largest sensible positive value for the scaler would 
be a value which would not cause overflow on the elements 
of C. 

The directive, X 

The directive, X, must be either an N, J, F, or L. When 

the directive is an odd number (J or L) and also d is 

less than 13, the rows of C will be punched in columnar 

form with a carriage return after each element. An N 

will be punched at the ent of each row and a J at the 

end of the matrix. 

When the directive is an even number (N or F) and also 

d is less than 13, the rows of matrix C will be printed 

across the teletype page. A carriage return and two delays 



will be punched after each n elements where n is the 
integer part of the quotient, 70/(d + l). At the end 
of each row of C, an N followed by two carriage returns 
and two delays will be punched. A J will be punched at 
the end of the matrix. 

Whenever d is 13 or greater, matrix C is stored on the 
drum, and no distinction is made between odd and even 
directives • 

Type T and Type D 

The directive at the end of a Type T parameter tape is 
either an N or a J. When a Type T tape is read, the 
computer will step on 2ko60. If the black switch is raised, 
a B matrix will be read from tape . (See Note 3) 
The directive at the end of a Type D parameter tape is 
either an F or an L. When a Type D tape is read, the 
computer will stop on 2^03 J. If the black switch is 
raised at stop 2^0SJ, a B matrix will be read from the 
drum . This matrix will be in fact the previously formed 
product matrix of size i by k stored at drum location D. 
The matrix will be read a row at a time, rescaled, and 
restored as a B matrix at location 2560 . The method of 
use for this is explained under the heading CASE III. 
It is also possible to specify the size and drum location 
of a matrix which is to be treated as a B matrix. After 
a Type D parameter tape has been read and the machine 
stops on 240SJ, if the white switch is raised up and down, 
three additional parameters will be read as follows s 

D space drum location of matrix 

r space number of rows 

c space number of columns 
To exercise this option, the matrix must have been stored 
on the drum at any time previously in successive locations 
by rows with a sum check at the end of each row. This is 
illustrated under the heading CASE IV. 
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ILLUSTRATIVE EXAMPLES OF PARAMETER TAPES s 

Type T Both B and A matrices are read from tape 
Parameters Meaning 

10 space N Print the elements of C to 10 decimal 

places by rows across the teletype page. 

k space +2J Multiply each element of B by 100 prior 
to the matrix multiplication; print the 
elements of C to k decimal places by 
rows with a carriage return after each 
element . 

k^8 space N Store C by rows with a row sum check at 
the end of each row in successive loca- 
tions starting at drum location k^8. 

Ty P e D The B matrix is read from the drum, A matrix from 
tape 

6 space -IF Read rows of previously formed pro- 
duct from the drum, multiply each element 
by l/lO, and store as a B matrix 3 form 
the new product C and print each element 
to 6 places J print by rows across the 
page. 

6^00 space F A matrix ML pi , is read from drum 
wh. switch location ^ D> 

6100 space 6100, and restored at 2560 as a B matrix; 
53 space matrix A is read from tape and the pro- 
21 space duct matrix C is stored by rows with a 
sum check at the end of each row in 
successive locations starting at 6k00. 

METHOD OF USE: CASE III, A (p) . ,.A (2) A (l) B = C st ops 

It is necessary first to form the product, A^ B=C^, 

and to store this on the drum at location D. 

1. Master tape 3402L 

2. Parameter tape, Type T 2^060 

3. Data tape B 2^0J7 
k. Data tape A^ ' (j term, sym.) 2^02L 
Next a Type D parameter tape is read. If it is not 
necessary to change the scaler or to change the store 
location on the drum, the parameter tape can be followed 
by several A matrices terminated by L symbols. 

5» Parameter tape, Type D 2^0SJ 

6a. Data tape A* ' (L term, sym.) 240SJ 
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6b. Data tape A 
etc. 



(3) 



(L term, sym.) 



240SJ 



METHOD OF USE! 



If it is necessary to change the scaler or the drum 
store location, or, if it is the end of the sequence, 
and it is desired to print the final product, an A 
matrix can be terminated with a J, which enables the 
user to read a new parameter tape. 
6x. Data tape k^' 1 ' (J term, sym.) 2^02L 

7. Parameter tape, Type D (with "d") 2*1-03 J 

8. Data tape A * p ' (J term, sym.) 2^02L 
CASE IV, A. r M rc = C. c Stp^s 
Let us assume that a series of steps (perhaps CASE I or 
CASE II ) formed the matrix, M , at drum location D* 
by rows with a sum check at the end of each row. Sub- 
sequent operations have cleared the Williams memory but 
hot the drum. To form the product, A ir M rc , the steps 
are as follows g 

1. Master tape 3^02L 

Parameter tape, Type D 240SJ 



/% r- 



2b . Parameters : 
, •* '■"" 
D ; space r space c space 

Raise white switch up and down 



5. Data tape A (with J term. sym. ) 



2ltf)S7 
2402 L 



STOPS AND ERROR DIAGNOSIS: 



Stops 


Loc. 


Diagnosis 


3^02L 


18J 


Master tape has read correctly^ bl. sw. 
to read par. tape. 


2^060 


053 


End Type T par. tape 5 bl. sw. to read B$ 
wh. sw. to convert to a Type D tape. 


2ifOSJ 


053 


End Type D par* tapej bl. sw. to restore 
B and read A| wh. sw. to read D*, r, c. 


21J-0SF 


05L 


D*, r, c, have been read; bl. sw. to re- 
store M and read A; wh. sw. to read B 
from tape* 


2^0JT 


08S 


End of tape Bj'-blo sw. to read A. 


2402L 


OJJ 


End of tape A with J term, sym.; bl. sw. 
to read par • 


2^088 


o6s 


F term. sym. end of 1st row of B$ bl. sw. 
to continue. 
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Stops 
20060 

200F5 
FF001 
FF002 



FF003 
FFOO^ 

FF005 
FF006 

FF007 
FF008 



Loc. 
087 

OJL 
06K 

oQk 



of8 

0S2 

10K 
109 
18F 
OJO 



STORAGE OF MATRIX Bs 



Diagnosis 



IMMMMMaMW 



F term. sym. end of subsequent rows of Bj bl. sw. to 
continue . 

F term. sym. end of row of A| bl. sw. to continue. 

J or L end of 1st row of B5 wh. sw. to continue. 

Number of elements in rows of B are not equals wh. 
sw. to continue as if =. 

Row vector of A does not conform to column vector of 
B5 wh. sw. to new par* tape. 

Improper drum address for store of C5 wh. sw. will 
calculate first available address and store C. (See 
Storage of C on Drum) 

Overflew on element of C5 wh. sw. to read par u 

B sum check failure^ wh. sw. to try again. 

Master tape sum check failure; Reread. 

Row of C read incorrectly from drum; wh. sw. to try 
again „ 



In order that the columns of B can be read from the 
drum in minimum access time, the first element in 
each row is stored at 2^60 + (j - l)l„ I is determined 
by the size of k. The capacity for B using the entire 
drum is shown in the chart belowg 



STORAGE OF MATRIX C: 



k 


I 


max. ,j 


"< k-< 65 


65 


157 


u? < k < 129 


129 


79 


129 < k < 193 


193 


53 


193 < k < 237 


257 


39 


257 < k < 311 


321 


31 



If matrix C is net stored on the drum, i, the rows 
of A, is not limited, k and j can be determined from 
the chart „ In any event, the maximum k is 311 . 
The first available location for the storage of matrix 
C on the drum will be after the final element of matrix 
B. This can be determined from the following formula. 
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D = 256O + (j - l) I + k 
where I is the increment between rows of B. 
Matrix C is stored on the drum by rows with a row sum 
check at the end of each row. The number of locations 
required will be i (k + l) and 

D + i(k + l) < 12,800. 

To direct the computer to store C on the drum, the 
parameter, d, must be greater than 12, and d will be 
interpreted as the drum address, D, at which to store C. 
If d is any number greater than 12, but less than D, 
the first available drum address, the computer will 
consider this to be an error and will stop on FFOO^. 
If the white switch is raised up and down, the computer 
will print a 5-place decimal address which is D, the 
first available address . Matrix C will be stored in 
successive locations starting at D. 

For CASE III when multiple products are formed D must 
be chosen to be xurge enough to aecoranodate the largest 
B matrix. 
DURATION IN SECONDS ?, About 30 seconds for the reading of the master and 

parameter tapes • 

Time : depends upon options exercised which are listed 
below; ' - ;v ' 

la. To read matrix B from tape and store on drum 
jk (.005 d^ + ;0l3) 

lb. To read matrix C from drum and restore as a B 

matrix 

i (.005 k + 0OI8) 

2. To read matrix A and caleula+e 

ij (.005 d + «010) + ik {^COk j + .009 ) 
si 

3a . To punch matrix C 

ik ( 9 0l67 d + .03*0 + o0b7 i 

3b. To store matrix C on the drum 
i (.002 k + .008) 
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NOTE lj 



NOTE 2% 



NOTE 3°o 



Whenever matrix B is rescaled by a positive power of 
10, + S, the accuracy of C is reduced from 11 decimal 
places to (11 - S) places «, 

It is not necessary to reproduce a large tape to punch 
a J or an L terminating symbol. The final row 
terminating symbol N can be hand punched to an F. 
After the computer stops, insert a single J or L 
character in the reader and raise the black switch. 
A Type T parameter tape can be converted to a Type D 
tape merely by raising the white switch up and down 
at the stop 2^060. 



DAT E September 28, 19*39 
PREPARED B Yr^fthjjj "" 
APPROVED BY 




ns 



LOCATION 


ORDER 


NOTES PAGE 1 M 28 


Abs. Rel. S»m 


"1 












003K 


3 






OOF 00397F 




Location of .A vector 


k 






OOF 0O554F 




Location of B vector 


5 






OOF 00711F 




Location of C vector 


6 






OOF OOF 






7 






1«)F L521(H12) 
0017K 






17 




(A) 
(B) 
(O 


00S3 O0S3 
OOSlf OOS^ 
00S5 00S5 










00 


40S^ L52l(m2) 








(ML) 


80F OOF 


•by ll,3fc(£) 


1 

Counter for number of elements ,t 






(-) 


OOF 001000 0000 OOOOJ 

* 








(+) 


OOF 006250 000 


OOOOJ 


' - - 






(o) 


OOF OOF 










(1) 


00F001F 










(i-D 


001F 001F 










(10) 


OOF 0010F 








(12) 


OOF 0012F 










m 


OOF 0064F 










(70) 


OOF 0070F J 








(01) 


8511F OOF 










(IE) 


8611F OOF 


.... 








(D5) 


OOF 002560F 


. . 








0*} 

(fts) 


OOF OOF 
OOF OOF 


by k,lk 

by ^8 (P), 7(1$; 


Store of matrix C on drum 
) 






(06) 


OIF OOF 


1 








(Tl) 


L5F 40F 


by if, 8(B2) 








(12 )[ J015F 75F .. . 
(T3) N09F L5F 


Iby 5.(B2) 








(*) 


756F ibF 


p by 6(B2) 








(T5) 


NOF kOF 


by .5(B5) ; . 








(T6) 


S68(A2) L5F 


by 4(B5) , 








(T7) 


JOF jkF 


by 3, 5(B5) 


-■ ■"-"":,.' ■ '.:'.' 


— L 


1 

— ,_ 


W 


75F 40F [£ by 2(B2) 


■ " - \ — — 



LOCATION 



AIDS. 



^7 



67 



Rel. 




1 
2 

3 
k 

5 
6 

7 
8 

9 

10 
11 
12 

13 
ll* 

15 
16 

17 
18 

19 
20 
21 
22 

23 
2k 

25 
26 

27 
28 
29 



Sym, 



(T9) 
(T10) 

(P) 



ORDER 



801F L5F 
Nll^F L5F 
OOK 

*a^F klF 
8l4F 50F 
7^(10) S5F 
kOF 91^F 
321L kl(£k) 
L5F ^03F 
F0(12) 3213L 
F53F 50 (0) 
007F 40F 
L5(70) 66F 
'S5F 1032F 
42 (NL) L53F 
0020F 1*63(A4) 
2215L L53F 
IA(D2) 4o(d4) 
52^F 5015L 
26(N12) 405F 
L3^F 3230L 
L5^F 322 3L 
402F L5(-) 
4013F F52F 
402F 3230L 
5013F 7J(-) 
2620L 50(0) 
L5^F 0022F 
46l2(B4) Ll^F 
1+02F L5(+) 
4013F F52F 

402F 3230L 

5013F 7J(+) 
262 7L 4ll6F 



£by 2(B2) 
by 7(B2) 



NOTES 



PAGE 2 



M 28 



Read parameters 



from sum check 
6(A1) 



Store 
tests 



"d" at loc. 3 

d < 13 v 



Calculate N, no* elements per line 



} 



read scaler and directive jj 
Store at loc. 4 and 5 



Negative scaling multiplier at loc .13 



positive scaling multiplier at loc. 13 



LOCATION 


ORDER 


NOTES PAGE 5 M 28 


Abs, 


Rel. 


Sym. 


101F S3F 








32 






33 




3234L L5(l) 




test: directive is odd or even 




3^ 




42 (NL) L55F 


from 12,l4(Al] 


1 




35 




F0(l) 3236L 




test: directive is F or L 




36 




24 (Bl) 24 (B6) 




stop: 24060 - 240SJ 




37 




192F 401F 








38 




L5(10) 4243L 




white switch to read 




39 




4lF 8l4F 




D*, r, c. 


87 


40 
41 




50F 7^(10) 
S5F 40F 








42 




914F 3640L 


■ 


' 




^3 




L5F 40F 








44 




F543L 4243L 








45 




L51F L4lF 








46 




4oif 3639L 






s 


47 




L510F L4(D2) 


. 




1 


48 




40 (D5) 24l(B6) 


Stop: 240SF 








OOK 






96 





(Bl) 


50S4 50L 


from 36 (P) 


Read first row of B 




1 




26(N12) 401F 








2 




L52l(N12) 407F 

■ ■ 


N12 store order at loc. 7 




3 




LO(N) 1020F 








4 




4012F 504L 


.j 


V'.at loc* 12 - 




5 


' 


26 (B2) L5(l) 




tally "j " at loc . 11 




6 




4011F 4l6F 


\ 






7 




L51F L0(l) 


; • 


! 


8 




369L 2210 (B3) 


!' 




9 


' 


L0(l) 3210L 


■ ■ ' ' i 

IFF1: J or L at end first row of B 




10 




FF1F LO(l) 


fwh. sw. to continue 




11 




3610L 24lO(B3) 


Stop on F2. 24088 








OOK 


\ 


108 





(B2) 


K5F 4217L 


from 5(Bl) | Set test constants 


1 


1 




L5(C) L412F 


- 2(B6).': " I 


i I 5 




42 (T8) 42 (T9) 
427 (A3) L5(B) 


i . jj 

i 

* i 

i 4 


1 h 


! 


L412F 42 (Tl) 


! ■ 



LOCATION 



Abs 



tr 



126 



Rel. 



140 



145 



5 

6 

7 
8 

9 

10 
11 
12 

13 

14 

15 
16 

17 

o 

l 

2 
3 

5 

6 

7 
8 

9 

10 

11 

12 

13 


I 1 

2 

3 

5 
I 6 

7 



Sym. 



(B3) 



ORDER 



NOTES 



W 



P/tfffi 4 



,M,.,2fl,. 



ally 



42(T2) 42(T3) 

42 (T4) 42 (T5) 

42(T10) 0020F 

46 (Tl) L512F 

L0(l) 402F 

F5(64) 408F 

L52F FO(64) 

3613L 2215L 

L4(l) 402F 

F58F L4(64) 

2210L L5(D2) 

L4(D3) L08F 

409F 22F 

00K 

50S 4 ^OL 

26(N12) L0(l) 

367L F511F 

4211F L52l(N12) 

L07F 40F 

L3F 3210L 

FF2F 2210L 

L0(l) 328L 

2212L LO(l) 

3212L 202L 

50(0) 5010L 

26 (B4) 22 (B3) 

50(0) 5012L 

26(B5) 24 (Al) 

OOK 

K5F 4225L jfrom ll(B3),2l(B6) 

L5(b) 42 7L 

467L 4211L 

4212L 4217L 

4222L L34F f 

l 

3P15L L54F I 

36UL 5013F \ 

75F 40F \ 



Calculate drum increment, 1, 
for store of B at loc • 8 



preset drum store order at loc. 9 



Read subsequent rows of B 



"j" at loc. 11 



tests elements of rows are equal? 
FF2i elements /=; wh. sw. to 
continue 

Stop on F: 20080 



from 8,11, (Bl) 



Stop on J or L: 240J7 



Rescale rev of B and 



store on drum 



negative res ca ling 



Abg. 



LOCATION 
Rel. 

8 

9 
10 
11 
12 

13 
Ik 

15 
16 

IT 
18 

19 
20 
21 
22 

23 
2k 

25 



166 





1 

2 
3 

5 
6 

7 
8 

9 

10 
11 
12 

13 
14 

15 



JJvm. 



(B5) 



ORDER 

L57L 1^(1-1) 
407L LO(Tl) 
326L 2215L 
5013F T5F 
OOF 40F 
F511L 4211L 
4214L L0(T2) 
3611L L59F 
IA8F 4018L 
409F L5F 
OOF OOF 
F518L 4018L 
F517L 421?L 
L0(T3) 3217L 
L56F IAF 
k06F F522L 
4222L L0(t4) 
3622L 22F 
OOK 

K5F 4221L 
L5(A) iAllF 
0020F 467F 
k6 (T7) L5(B) 
iAllF 42 (T6) 
42(T7) 42(T5) 
4110F L5(D4) 
40(D5) L3(D4) 
3219L L316F 
3610L 2221L 
L5(D4) 409(A3) 
L0l8(B4) 3221L 
FF4F 92131F 
L5l8(B4) 409(A3) 
40(D5) L0(D2) 
40F 5022L 



NOTES 



by 25 (P) 



by 3L 
by 16L 
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M28 



positive rescaling 



Store on drum 



sum check for B at loc. 6 



from 13(B3),25(B6) 

Set additional test constants 



clear row counter for A at loc .10 



FF^-s Drum address too small 

for store of C; wh. sw. to 

find first available address and to 

Print corrected address 



LOCATION 



Abs. 



186 



189 



Bel. 

16 

IT 
18 

19 

20 

21 

22 



1 

2 

3 

k 

5' 

6 

7 
8 

9 

10 
11 
12 

13 
1^ 

15 
16 

17 
18 

19 

20 
21 
22 

23 
2k 

25 



Sym. 



+m 



(36) 



ORDER 



75F S5F 

505F 5017L 

26(P17) 92131F 

2221L 9259F froa 8L 

9259F 92139F 

92707F 22F 

OOF 00 0000 1000 0000J 

L510F tollF from 36 (p) 

kl6F 501L from k9{$) 

26 (B2) L5(D5) 

L0(D6) ^o8l 

10(b) lM2F 

4217L kUkF 

4115F L5(B) 

1^29L tel3L 

OOF 00F by 3L 

40F 1*0F 

F58L 408L 

F59L ^29L 

F0(T5) 368L 

4ll^F L5F by 7L 

iAl^+F *K)ll*F 

F513L ^213L 

LO(TIO) 32131 

L514F LOF by 4L 

kOF L3F 

3220L FF8F 

26lL 502OL 
26(Bif) F515F 
i^215F ID11F 
322*fcL 226l 
OOF 5024L 

26 (B5) 26 (ai; 



NOTES 
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M 28 



Replace B with 



Read row of C from drum 
with sum check 



Sum check for row at loc. Ik 



FF8s Row of C read from drumj 
incorrectly 5 wh. sw to try again 



LOCATION 


ORDER 1 NOTES PAGET M 28 


AbSo J] 


tel. 


Sym, 


00K 

503 ^OL J 


. 




215 





(Al) 


From 13(B3),25(B6) Read rows of A 




1 




26(H12) LO(l) : 


L3(A3), lMAfc) 






2 




363L 2214L l 








3 




L0(1) 367L 








k 




L3(D^) 325L J 
226L 9283^F 






5 




Print "J" 




6 




92135F 2k(?) 




Stops End of A| 2^02L 




■7 




LO(l) 328L 








8 




201^L L3(D^) 




Stop on Fs 200F5 




9 




3613L L55F 








10 




101F LO(l) 








11 




323^(P) L5(D1) 








12 




kOl6¥ 223^ (P) 








1 x 




92834F 92135F J 




X X XU U V 




Ik 




223^ (P) L521(N12) 
1 






15 




L07F i^)F 




16 




L3F 36 (A2) 






233 


17 



(A2) 


FF3F 26 (P) 

OOK 

kU.k¥ L5(D1) 


from 16 (Al) 


FF3§ row vector of A does not con- - 
form to Bj wh. sw. to new parameter. 




l 




lA(D3) L0(1) 








2 




if09F F510F 




tally "I" at loc. 10 




3 




1*010F ^1F 








> 




L5(C) k22JL 








5 


. 


L5(B) te9L 








6 
7 




4213L F59F 
k09F kO&L 








8 




OOF OOF 




Read column of B from drum 




9 




kO¥ kOF 




■ 




10 




L58L IA8F 




''.-■■ 




11 




408L F59L 




" . ■ ■ ■ > 




12 




1^29L L0(T5) 




■ ' ■ ; ; ' 




13 




368L L5F 








ik 




iAl^F ^Ol^F 








15 




F513L 4213L 


1 f orm sum of b. . at loco ik 



LOCATION 


ORDER 


NOTES PAGE 8 M 28 


Abs. 


Relo 


Synu 


L0(T6) 3213L 








16 






17 




L5(A) 4620L 








18 




L5(B) 1+220L 








19 




kZF S5F 






253 


20 
21 
22 

23 
2k 

25 




50F 7^F 
401F 3235L 
L5F 323^L 
IA1F 3233L 
liOF L520L 
iA(l-l) >020L 


by 17, 18L 


vector multiplication 




26 




L0(T7) 3219L 








27 

28 




L5F iK)F 

- 

F527L te27L 


by 4L 


store C.' 
ik 




29 




L0(t8) 365L 




.'.'■."■ 




30 




L56F LOl^F 








31 




40F L3F 


. 


■ . 




32 




36(A3) FF6F 




FF6s B sum check failure | 




33 




26L FF5F 




wh. sw. to try again 




3* 




26 (p) Lto 








35 




2624L L5F \ 




FF5 § Overflow on element 




36 




3637L 223^L 




of C| wh. sw. to read a 




37 




IAIF 362J& 


.■.-■- 


new parameter 


272 


38 



W 


2233L OOF 


■ 
from 32 (A2) 


■ 




1 




L5(0) 423L 




Store row of C on drum with 




2 




k2&L klF 




row sum check at end of each 




3 




OOIF L5F 


. 


row. 




if 




I^F 40F 




...■"'■■' 




5 




F53L>23L 








6 




L0(T9) 323L 




■- ■ 


279 


7 
8 




L5F 40F 
OOIF L5F 


by 3(B2) 




1 


9 




OOF OOF 


by 10, 13 (B5) 


Store on drum 




10 


: 


F59L ^0% 

— ; ■ • 




1 



LOCATION 


ORDER 


NOTES PAGE 9 M 28 


Abs. 


ReTT" 


Synio 


F58L 428L 








11 






12 




F0(T9) 328L 






286 


13 


1 
2 


W 


22 (Al) OOF 
L5?C) 422L 
lOl^F OOIF 
lfll5F L5F 


from 0(A3) 


Print row of C 




3 
k 




50F 503L 
26(P17) F52L 


by 12 (P) 






5 




h2ZL F51^F 








6 




1+21^F L012F 








7 




36UL F515F 




tests end of row? 




8 




4215F 10(NL) 


■ ' 






9 




522L 9213IF 


' 


tests end of n elements 




10 




92519P 262L 








11 




92770F 92135F 


..■'■■ 






12 




92519F 22 (Al) 






299 




(K12) 


OOK 




Input Routine (N12) 


338 




(P17) 


OOK (EL7) 
OOK 




Print Routine (P17) 


397 






L3F 3^(P) 






398 






FF7F 26 (P) I 


FF7s Read failure on 


399 






Fif3691F 580^78F 


master tape G 






- 


26L 261N 
















■ ■ 



